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Abstract 

In the present work we apply the atomic approach to the single impurity Anderson model (SIAM). 
A general formulation of this approach, that can be applied both to the impurity and to the lattice 
Anderson Hamiltonian, was developed in a previous work (larXiv:0903.0139"V l [cond-mat.str-el]). 

£ : 

The method starts from the cumulant expansion of the periodic Anderson model (PAM) , employing 
the hybridization as perturbation. The atomic Anderson limit is analytically solved and its sixteen 

7— I ■ 

■ eigenenergies and eigenstates are obtained. This atomic Anderson solution, which we call the 

: 

, (AAS), has all the fundamental excitations that generate the Kondo effect, and in the atomic 

'xf ■ 

approach is employed as a "seed " to generate the approximate solutions for finite U. The width of 
I ' the conduction band is reduced to zero in the AAS, and we choose its position so that the Friedel 

o\ : 

<^ . sum rule (FSR) be satisfied, close to the chemical potential /x. 

• i-H ! We perform a complete study of the density of states of the SIAM in all the relevant range 

of parameters: the empty dot, the intermediate valence (IV-regime) ,the Kondo and the magnetic 
regime. In the Kondo regime we obtain a density of states that characterizes well the structure of 
the Kondo peak. To shown the usefulness of the method we have calculated the conductance of a 
quantum dot, side coupled to a conduction band. 

PACS numbers: 36.20.Ey,64.60.Cn,05.50.+q 
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I. INTRODUCTION 

Experimental results indicate [2j that localized magnetic moments could appear when 
magnetic impurities are dissolved in a metal, and a minimum in the electrical resistivity 
is always present in the material at low temperatures T. Jun Kondo was the first who 
associated the formation of magnetic moments in the material with the minimum of the 
resistivity. He studied the phenomena employing perturbation theory in second order in 
the Born approximation, [3] and showed that the minimum was generated by the spin- 
flip scattering of the conduction electrons with the magnetic moment of the impurity; this 
phenomena become known as the Kondo effect. In its calculation Kondo showed that the 
resistivity increases with the decrease of temperature according to ln(T), which become 
known as the logarithmic divergence of the Kondo effect. In the alloys formed in such a way 
it is not possible to control the Kondo parameters microscopically, but with the advent of 
nanotechnology the Kondo effect was realized experimentally in quantum dots (QD), with 
complete control over all the relevant parameters Q]. 

Quantum dots are artificial atoms that offer a high degree of parameter control by means 
of simple gate electrostatics, allowing the nanometric confinement of the electrons. The 
coupling between the quantum dot and the electronic reservoirs produces tunneling events 
that can modify the discrete energy levels of the dot, changing them into complicated many- 
body wave functions. When conducting electrons move in an out the nanostructure, a 
progressive screening of the atomic spin occurs, in complete analogy with the well-known 
Kondo effect in solids containing magnetic impurities. The Kondo effect in quantum dots is 
then observed as a zero-bias conductance resonance, associated with the entangled state of 
the electrons in the electronic reservoirs and in the dot. 

The experimental realization of the Kondo effect in QD's renewed the interest in search- 
ing new approximate solutions for the single impurity Anderson model (SIAM) and several 
calculations of the Green's functions (GF) employing the method of motion equation (EOM) 
were reconsidered, jf], |(| but such methods suffered several drawbacks due to their failure to 
satisfy the completeness and the Friedel sum rule. 7J Another method that was extensively 
applied to describe the SIAM is the slave boson mean field theory, s| but this method 
suffered from a spurious second order phase transition. |9( Another recent approach that 
describes the SIAM is the local moment approach, pLOj] that satisfies a number of correct 
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limits: weak-coupling Fermi liquid, atomic limit, and the Kondo regime, but the strong 
coupling limit is plagued with a symmetry breaking and an unphysical static order value 
of the local magnetic moment. Very recently appeared an analytical calculation, 11( that 
employed diagrammatic expansions with Feynman diagrams and obtained the Kondo reso- 
nance in the electron-hole symmetric limit, but out of this limit the Kondo resonance was 
lost. Another appealing approach is the quantum Monte Carlo technique, [12j that is ex- 
pensive from the computational point of view and is restricted from intermediate to high 
temperatures. Other expensive computational numerical methods are the dynamical mean 
field theory (DMFT) 13] and the density-matrix renormalization group (DMRG). [ljj] The 



most effective approach which nowadays describes t 



the numerical renormalization group (NRG) 
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re dynamical properties of the SIAM is 
17| but as the other powerful methods 



it is expensive from the numerical point of view and presents some limitations; the NRG 
describes very well the Kondo Physics but presents convergence problems in the extreme 
Kondo limit. [3] Until now it was missing a dynamical theory that describes the transition 
from the weak coupling limit U « A to the strong coupling limit U » A, where U is 
the Coulomb repulsion and A = irV 2 p c , and V is the hybridization and p c is the density of 
states of the conduction electrons. 

The main objective of this paper is to fill the gap described above: the atomic approach 
is able to describe quantitatively all the regimes of the SIAM in the weak, intermediate 
and strong correlation limits of the model. Due to the simplicity of its implementation 
(practically all the method is analytical) and very low computational cost (a density of states 
curve can be obtained in few seconds or less), the atomic approach is a good candidate to 
describe strongly correlated impurity systems that exhibits the Kondo effect, like quantum 
dots [19( or carbon nanotubes. 20] The atomic approach does not substitute powerful 
computational methods like NRG, DMFT or DMRG, but can be a first choice to describe 
systems where the Kondo physics is relevant. In an earlier work 2l|] we developed the atomic 
approach for the strong coupling limit (U — ► 00), and in the present work we present the 
atomic approach for finite U. 

The Kondo effect for both the impurity and the lattice are described by the Anderson 
Hamiltonian, which considers systems that have two types of electrons: the localized elec- 
trons (/-electrons), that are strongly correlated, and the conduction electrons (c-electrons), 
that can be described as free. The Anderson model allows the interchange between the / 
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and the c electrons through the hybridization interaction, and to study this problem we 
shall consider the cumulant expansion of the periodic Anderson model (PAM) employing 
the hybridization as perturbation. 22| From the cumulant expansion we obtain formal 
expressions for the exact one-electron Green's functions (GF) in terms of effective cumu- 
lants, that are as difficult to calculate as the exact GF, and our approximation consists in 
substituting these effective cumulants by those of the atomic case of the Anderson model, 
that is exactly soluble. In general the Anderson model does not have analytical solutions, 
but when the energies of all the N conduction states have collapsed (the conduction band 
has zero width) and the hybridization is local (i.e. independent of the wave vector k), the 
Anderson Hamiltonian has an exact solution, and the electronic Green's functions can be 
calculated analytically. In this case it is only necessary to study a single site, that can 
hold up to four electrons: two /-electrons and two c-electrons. With the /-electrons we 
can build four states: one | 0) with no electrons, two with one electron with spin up| +) or 
with spin down| — ) and one | d) with the two electrons at the same site. In a similar way 
we can build up four states | 0), |f), |TI) with the c-electrons, and for a given site we 
have then a vector space of dimension sixteen. The corresponding Anderson Hamiltonian 
is then easily diagonalized, giving sixteen eigenenergies and eigenstates that we shall name 
the atomic Anderson solution (AAS), and this minimal Anderson system contains already 
the fundamental excitations that generate the Kondo effect. To chose the position of the 
conduction band of the AAS, we impose the satisfaction of Friedel's sum rule. 

In Sec. [Til we present a brief review of the basic equations of the atomic approach formal- 
ism developed in an earlier publication. [l| In Sec. [Ill] we discuss in detail the equations and 
approximations employed in the development of the atomic approach. In Sec. IIVI we define 
the completeness and calculate the occupation numbers. In Sec. |V]we discuss what criteria 
we shall use to determine the parameters of the AAS: the satisfaction of the completeness 
or the Friedel sum rule. We also present a discussion of the emergence of the Kondo peak 
from the weak to the strong coupling regime. In Sec. IVII we present the results of the 
formation of the Kondo peak for a fixed Ef level and variable correlation U. In Sec. IVII I 
we fix the correlation U and vary the Ef level describing the empty dot, the intermediate 
valence (IV), the Kondo and the magnetic regimes. To show the usefulness of the atomic 



approach, we calculate in Sec. IVIIII the conductance of a side coupled QD. [23|, |24| In Sec. 
I Villi we present the conclusions of the work and in the appendix |A] we give the details of the 
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analytic calculation of the exact Green's function in the limit of zero conduction bandwidth, 
as well as the relevant atomic Green's functions. 



II. THE ANDERSON HAMILTONIAN 

In this section we present a brief review of the basic equations of the atomic approach 
formalism developed in an earlier publication The Hamiltonian for the Anderson lattice 
with finite U is given by 

+ U Y,j n j,v n j,a + H h , (1) 

where the operators C£ a and Ck jCr are the creation and destruction operators of the conduc- 
tion band electrons (c-electrons) with wave vector k, component of spin a and energies E^ a . 
The fl a and /j jCr are the corresponding operators for the /-electrons in the Wannier localized 
state at site j, with site independent energy E a and spin component a. The third term is 
the Coulomb repulsion between the localized electrons at each site, where rij tU = f] a fj^ is 
the number of /-electrons with spin component a at site j and the symbol a denotes the 
spin component opposite to o. The fourth term describes the hybridization between the 
localized and conduction electrons 

H h = }.,,,/:/',,,, + v;. k .Jt, f,,,). (2) 

The hybridization constant Vj^ a in this equation is given by 

^W = ^^(k)exp(2k.R i ), (3) 

and when the Hubbard operators are introduced into Eq. ([2]) , the hybridization Hamiltonian 
Hh is transformed into: 

H h = ^ \Vj ba>ka Xj ba C k(T + V* baMa Cl a X jM ^j , (4) 

j6a,k(T 

where the label a = (b,a) in describes the transition | a >— >| b >, and the local state 
| a > has one electron more than the state | b >. There are four local states | >, | + >, 
| — > and | d >=| +, — > per site, and there are only four X operators that destroy one 
local electron at a given site. 



The identity relation in the reduced space of the localized states at site j is 



Xj, o + Xj )aa + + Xj iC id — I, (5) 

where a = —a, and the four Xf m are the projectors into the corresponding states | /, a). The 
occupation numbers on the impurity n/ >a =< X^ aa > satisfy the "completeness" relation 

ra/,o + n f,a + nf,w + rif,d = 1- (6) 
We use the index I x = 1, 2, 3, 4, defined in Table I, to characterize these X operators: 
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1 


2 


3 


4 


a = (b, a) 


(0,+) 


(o,-) 




(+,d) 



TABLE I: Representation of the possible transitions present in the finite U atomic SIAM Hamil- 
tonian. I x = 1,3 destroy one electron with spin up and I x = 2, 4 destroy one electron with spin 
down. We use a = + and a = — instead of a =T and a = j to emphasize that the spin belongs to 
a local electron. 



To simplify the calculation we now introduce the two matrices 

{M} aa ,=M e JJ(k,z,u), (7) 

and 

{W} , =Wa>, a (k,a,z), (8) 
r-| 

where M is the effective cumulant matrix |l( an d the matrix elements of W employed in 
the PAM calculation are defined by 

W^ a (k,a,z n )=V(a',k,a)V*(a,k,a) g% (k, z n ) , 

where z n = iu n are the Matsubara frequencies and 

is the free GF of the conduction electrons. A related matrix appears in the impurity case 

{W} ; =Wa>, a (<T,z). 



with matrix elements defined by 

W a ,, a (a, z) = ±r V(a/, k, a)V*(a, k, a) G% (k, z) . 
s k 

The hybridization is spin independent in the Anderson model, so then 



V(0<x, k, a) = V(ad, k, a) = 

V(0a, k, a) = V(ad, k, a) = 0. (10) 

We shall assume a mixing that is only local, so that V^-(k) in Eq. is k independent, 
and in Eq.(jH) we then have 

V(0a,k,a) = V, (11) 

V(ad,k,a) = aV, (12) 

where we have also assumed that V^-(k) is independent of a — ±1. 
We shall use below the matrix A = W.M, with matrix elements 

A aa , (k,a,z) = (W-M) aa ,= 

W aai (k, a, z) M aia , (k, a, z) , (13) 

and when the Hamiltonian is spin independent or commutes with the z component of the 
spin, the 4x4 matrices G** , M, W and A can be diagonalized into two 2x2 matrices, 
e.g.: 

In this matrix the indexes I x defined in Table I have been rearranged, so that I x = 1,3 
appear in G* and I x = 2,4 appear in G, . 
Employing Eqs. ffTTlfT2l) we find for the PAM 




w T (M = \v\ 2 g° cA (k,z) | " - 1 , (is) 

i -i 



Wi(M = \v\ A g° Ctl (k, z ) | i i ] , (i6) 



where Q® a (k, z) is given by Eq.Q. For an impurity located at the origin we find instead 



W T (s) = \V\ 2 ^{z) ( | | | , (17) 



•1-1 



W^z) = \V\ z n (z) ^ i | , (18) 
where 

^) = 4-E C(M- (1 9 ) 

s k 

For a rectangular band with half-width D in the interval [A, B], with B = A + 2D we then 
find 

K.M^fcfc^Y (20) 



2D V^ + ^ + Z 1 , 

where the chemical potential /i appears in y>o-(z) because of the e{k,a) = E k a — \i in 

SS„(M.. 

From our earlier work |l| we obtain the exact formal Green's functions, both for the PAM 
and for the SI AM: 

G^M^I-A*)- 1 , (21) 

and from this equation follows 

M CT =(I + G^-W CT ) _1 -G^. (22) 

III. THE ATOMIC APPROACH FOR THE SINGLE IMPURITY ANDERSON 
MODEL (SIAM) 

Defining the exact cumulants as 

M t =("A ; M, J"*), (23) 
\m 3 i m 33 J \m 4 2 m 44 J 

one obtains the exact Green's functions Gy(iu) of the localized / electrons by performing 
the matrix inversion in Eq. ( l2~Ti) . employing Eqs. ( TT71fT8l) and Eq. ( 1231) : 
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G{'(iw) 



G{'( 



IU) 



rff rff 
r ff rff 

U 31 ^33 



rff rff 

^22 ^24 

rff rff 

^42 ^44 



m 31 m33 



\V\ 2 <fi{(iuj)) (m u m 33 - m 13 m 31 




m 22 m 2 4 

m 42 m44 



\V\ ip^((iu})) (m u + m 33 + m 13 + m 31 ) 



\V\ 2 (fi((iuj)) (m 22 m 44 - m 24 m 42 




1 _ 1^1 fldiv)) (m 22 + m 44 - m 24 - m 42 ) 



(24) 



■ (25) 



In the same way we can obtain the conduction G£ c ( 
Green's functions; a detailed derivation can be found in 



i, k',iu) and the cross G^(k, iu) 
1 



\V\ 



Gf (k, k', iu) = G° cA (k, iu) 5 (k, k') + 

(mn + m 33 + m 13 + m 31 ) . 

2 (k , za;) , 

1 - |V| (^ T (iu;) (m n + m 33 + m 13 + m 3 i) ' 



(26) 



I VI 



<i (k, iw) 



Gf (k, k', iu) = g° cd (k, iu) 5 (k, k') + 
(mm 22 + m 44 - m 24 - m 42 ) 



-<??, (k',iu,), 



(27) 



N s C,ly ' 1 - | | 2 ^l(^) ( m 22 + ™44 - "^24 - "^42) 

and the cross Green function G C J is defined by a column vector with two components as [l| 



G| (k, iu) 



Gf (k, iu) 



G C J 



r c f 
r°f 



m u + max , ™i3 + ^33 



V 



|V| ¥? T (io;) (mn + m 33 + m 13 + m 3 i) 
m 22 - m 42 , m 24 - m 44 



(28) 

(29) 
(30) 



1 y v^Vl C ' i V 7 1 - \V\ 2 tp[(iu) (m 22 + m 44 - m 24 - m 42 ) 

The calculation of the exact effective cumulants M a is as difficult as that of the exact 

Gf/ , and the atomic approach consists in using instead the effective cumulants of a similar 

model that is exactly soluble. The atomic limit of the SIAM is just a particular case of 

the general model, and we shall then use the AAS to calculate the exact Green's function 

GfJ' at (z) of the atomic problem, which then satisfies a relation of the same form of Eq. ( l2~li) : 
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From this equation we obtain the exact atomic cumulant M"* 

Mf= (I + G f J' at ■ W°) ^ ■ G f J' at , (32) 
which satisfies Eq. (1221) . For an impurity located at the origin we change Eqs. (11 711 181) into 

W°(z) = |A|V(z) r ij , (33) 
Wl(z) = \A\ 2 rt(z) ^ , (34) 

where 

CO = _1 , (35) 
z — e Q — /i 

is obtained by replacing all the Q® a (k, z) in Eq. (fT9l) by those corresponding to the zeroth- 
width band located at Eq, namely the bare conduction Green function. This procedure 
overestimates the contribution of the c electrons, 29| because we concentrate them at a 
single energy level e Q , and to moderate this effect we replace V 2 by A 2 in Eqs. (I33H34I) . 
where A = irV 2 /2D is the Anderson parameter. 

The atomic approach consists in substituting the exact effective cumulant M a , that ap- 
pears in Eqs. ( |24ll30|) . by the exact atomic one M**, which is defined by Eqs. ( I32ll35l) . We 
call M" p this approximate cumulant, and we make the substitution 

M a -> M7 (36) 

Performing the matrix inversion in Eq. (1321) and employing Eqs. ( 13311341) it is now straight- 
forward to obtain 



m ll m 13 i V 931 ^ 33 



911 913 I , I A|2 of W 

+ |A| ip° } {iLo) {giig 33 - 913931 




M TM = 1 1 = ~ i L -TT^7— l — : : — ^ (37) 



an ap 
m 31 m 33 



1 + | A| (gn + #33 + 9i3 + 93i) 



922 924] 2 (11 

+ |A| (flilU) (922944 - 924942) 

, ™Z ™24 \ \ 9 ^ 9 ^l V 1 l , 

M* p (iu) = 22 24 = -A '- , ^ (38) 

mZ mZ 1 + I A l Vliiu) (922 + 944 - 924 - 942) 
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where the g^- are the atomic Green's functions of the / electrons that are calculated in 
Appendix El Substituting now the approximate M" p 's in Eqs. f l24ll25|) we obtain the Green's 
functions of the localized / electrons in the atomic approach. 

We should stress that it is essential to use ^° a (z) rather than (p a (z) in Eqs. (I33|I34I) to 
obtain a well defined Kondo peak structure at the chemical potential /i. If we perform instead 
the calculation employing (p a (z) we always obtain a wrong structure with two or more peaks 



around the chemica 
Anderson model 27 



potential, as obtained in early works using the atomic solution of the 



281 ] . In the computational calculation we fixed the chemical potential 



at fi = and varied the conduction atomic level Eq in such a way that the Friedel sum rule 
should be satisfied. This point will be discussed in more detail in the next section. 



IV. THE COMPLETENESS PROBLEM AND OCCUPATION NUMBERS 

Employing Eq. ( fl4l) we separate the different quantities in two different types: those 
associated with I x = 1,3 corresponds spin up electrons and those associated with I x — 2, 4 
corresponds to spin down electrons. 

In the finite U case, the identity operator in the space of the impurity local / states is 
given by 

X OQ + X aa + X-^ + Xdd = I- (39) 
The completeness is associated to the average of this equation: 

(X 00 + X aa + X^ + X dd ) = 1, (40) 

where the first term is the vacuum occupation number, the second and the third terms 
are the occupation of the spin up and down respectively, and the last term is the double 
occupation. Employing the notation in Table I to identify the X^ a operators we could then 
write this equation in the form 

(X^xt) + (xlXj) + (x 3 Xl) + (X*X 3 ) = 1, (41) 

and all the different averages could be calculated employing the Green's functions G{[{uo) 
and G{1 (uj) defined in Eqs. (1241) and associated with the processes I x = 1, 3: 

< X 0<0 >=< X X X\ >= (^-j J duIm(G{{){l - n F ), (42) 
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< X +i+ >=< X\X X >= (^pj J dcuIm{G{^)n F , (43) 
< >=< X 3 Xj >= (^pj J dwIm{G(i){l - n F ), (44) 

< X d4 >=< XlX 3 >= f duIm(G f 3 ()n F , (45) 

where n F (x) = 1/ [1 + exp(j3x)] is the Fermi-Dirac distribution. 

In a similar way we could employ the Green's functions G 2 2 an d G{[ defined in Eqs. 
(125]) . and associated with the processes I x = 2,4: 



X 2 X\) + (X 2 f X 2 ) + (x 4 X 4 f ) + (x 4 f X 4 ) = 1, (46) 

< X . >=< X 2 X\ >= (^j J duIm(GU)(l - n F ), (47) 

< X_,_ >=< XjX 2 >= (^pj j duIm(G f i)n Fl (48) 

< X +i+ >=< X 4 X] >= f^J J dujIm{G{{){\ - n F ), (49) 

< X d4 >=< X|X 4 >= (^pj j duIm(G{()n F . (50) 

V. THE FRIEDEL SUM RULE: A CRITERIA TO BE SATISFIED 



The Friedel's sum rule (FSR) |30j gives at T = 0, a relationship between the extra states 



induced below the Fermi level by a scattering center and the phase shift at the chemical 
potential i] a (fi), obtained by the transference matrix T// ;(J (z) = V 2 G l J^ p (T (z), where V is the 
scattering potential. For the SIAM the extra states induced are given by the occupation 
number n/ jCr of the localized state, and the scattering potential is the hybridization that 
affects the conduction electrons. The Friedel's sum rule (FSR) for the Anderson impurity 



model can be written as 



3l| 



sin 2 (ixn f(7 ) 

A7T 



where p/ jCr (/z) is the density of states of the localized level at the chemical potential. 

The atomic approach consists in substituting the exact effective cumulant M CT , that ap- 
pears in Eqs. ( I2"4"ll30~i) . by the exact atomic one M"*, which is defined by Eqs. ( I32ll35|) . The 
difference between the exact and the approximate GF's is that different energies Ek appear 
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in the c-electron propagators of the effective cumulant M^*(z), while these energies are all 
equal to the atomic conduction level Eq in M^(z). Although M^(z) is for that reason only 
an approximation, it contains all the cumulant diagrams that should be present, and one 
would expect that the corresponding GF would have fairly realistic features. One still has 
to decide what value of Eq should be taken. As the most important region of the conduction 
electrons is the chemical potential /i, we could choose £o = /i, but we shall consider instead 
that the position of the atomic conduction level is given by £ = fi ± 5eq . This leaves the 
freedom of choosing e so that the Friedel sum rule given by Eq. (15T!) should be satisfied. 



In our earlier paper 



21] , where we developed the atomic approach of the Anderson 



impurity model for infinite U, we imposed the satisfaction of the completeness relation 
Xj tOQ + Xj t(T(T + Xj tWW = I instead of the fulfillment of the Friedel sum rule. To be rigorous 
the FSR is only valid at temperature T = 0, but we can use it as an approximation at 
temperatures of order of the Kondo temperature T#. The validity of the completeness 
condition is much more general than the FSR, because completeness is valid for the whole 
range of temperatures and parameters of the model. In the case of the finite U Anderson 
model the completeness is always satisfied and we cannot use it to determine Eq: we shall 
then use the fulfillment of the Friedel sum rule as a condition to obtain the adequate physical 
solution. 

Employing the results of the Section [Til and Appendix [X] we can calculate the (/, c) 
components of the matrix Green's functions for finite U in the atomic approach, and the 
corresponding spectral densities are given by 



\-^\ Im[G ff > cc (Lu)]. (52) 



In all the figures of the paper we employ A units, with A = nV 2 /2D = 0.01, with 
D = 1.0. 

In Fig. [I] we calculate the density of states at the chemical potential /i as a function of Ef, 
for U — > oo and T = 0.001A. In this case we impose the satisfaction of the completeness 
and the figure shows that in the extreme Kondo region n/ ~ 1, there is a considerable 
departure of the pf (/i) that satisfies the Friedel sum rule; in the inset of that figure we show 
the evolution of the occupation n/ toward the Kondo limit and that the completeness is 
satisfied. 

In Fig. [2] we impose the satisfaction of Friedel'sum rule. According to that figure the / 
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FSR 




nl i I i I i I i I i I i I i I i I i I i I 

-18 -16 -14 -12 -10 -S -6 -4-2 2 

E/A 

FIG. 1: Density of states at the chemical potential p as function of Ef, for U — » oo and T = 0.001A; 
we impose the satisfaction of the completeness relation. The pf (p) in FSR was calculated with 
Eq. HEED- 
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E/A 

FIG. 2: Density of states at the chemical potential /x as a function of Ef for U — > oo and T = 
0.001 A. The pf (p) in FSR was calculated with Eq. (15 II) . In this case we impose the satisfaction 
of Friedel'sum rule. 

density of states at the chemical potential /z, satisfies the FSR. In the inset of this figure we 
show that the completeness is lost but in the Kondo region the departure from completeness 
is very low; the error is less than the 1%, which justifies the use of FSR to determine Eq. We 
shall then employ this criteria both for U — > oo and for finite U, because the temperatures 
involved in this effect are generally very low. 

Next we study the evolution of the Kondo peak for finite U when the correlation energy U 
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U=50.0A |1 

U=100.0A ii 
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U=Infinite > ' 



o 
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FIG. 3: Density of states for T = 0.001 A, E f = -10.0A and several U values: U = 20.0A, 
U = 50.0A, U = 100. OA, U = 500. OA and U = oo. 
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FIG. 4: Detail of the lower band in the density of states of the Fig. El 



increases toward U — » oo. In Fig. [3]we plot the density of states corresponding the evolution 
of the Kondo peak for T = 0.001A, E f = -10.0A and several U values: U = 20.0A, 
U = 50. OA, U = 100. OA, U = 500. OA and U = oo. We show the upper band only for 
U = 20. OA. In Fig. H] we plot in detail the resonant band located around Ef: the curve 
with U = 500. OA is practically coincident with the one with U — > oo. In Fig. [5] we plot in 
detail the Kondo peak for the same values, and again the plot with U = 500. OA practically 
coincides with the one for U — > oo. 
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FIG. 5: Detail of the Kondo peak in the density of states of the Fig. 



VI. THE EMERGENCE OF THE KONDO PEAK 
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FIG. 6: Density of states p f (u) for U = 1.0A, E f = -10. OA and T = 0.001A. 

In the set of Figs. ([B UTT]) we show the evolution of the / density of states as a function of 
the Coulomb repulsion U. We begin with the uncorrelated limit of the finite U SIAM (weak 
coupling limit), and we consider U = 1.0 A, U = 5.0 A, U = 10.0 A, U = 15.0 A, U = 20.0 A 
(this value corresponds to symmetric case) and U = 25. OA. 

We observe that in the Figs. ([MED, where the correlation is weak, the three peak structure 
characteristic of the SIAM starts to appear, but in this region the Kondo peak is not yet 
formed. However, for U ~ 15. OA we are already in the Kondo regime for finite U, as shown 
in Fig. [91 where the Kondo peak is well defined. The Fig. [10] for U = 20. OA corresponds 
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FIG. 7: Density of states p f (u) for U = 5.0A, E f = -10. OA and T = 0.001 A. 
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FIG. 8: Density of states p f (uj) for U = 10. OA, E f = -10.0A and T = 0.001A. 

to the symmetrical limit of the Anderson model, in which the total occupation number is 
exactly n/ = 1.0. It is interesting to observe that in this case the atomic approach is able to 
reproduce the correct symmetry of the density of states. Finally in Fig. [II], for U = 25. OA 
the Kondo peak continues pinned to the chemical potential. 

In Fig. [12] we present the / density of the states at the chemical potential /x. At the 
uncorrelated side (U < 15. OA), the density of states at the chemical potential is small, 
but increases as the value of the correlation U increases. For (U > 15. OA) we attain the 
Kondo regime, where the occupation number nf ~ 1 and the Friedel sum rule produces 
(p/(/i) = l/VA) |32j. In the inset of the figure we present the density of states of conduction 
electrons c at the chemical potential /i, showing a loss of conduction electron states that 
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FIG. 9: Density of states p f (uj) for U = 15. OA, E f = -10.0A and T = 0.001A. 
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FIG. 10: Density of states for T = 0.001 A , U = 20A and E f = -10.0A. This parameter set 
corresponds to the symmetric case of the model and the total localized occupation number is 
exactly the unity rif = 1.0 

migrate to the localized band, screening the impurity and generating the Kondo effect. 

In Fig. [13] we present the total localized occupation number rif — n / )Cr as a function 
of the correlation U in A units. At the non correlated side (U < 15. OA), the / occupation 
number assume values between rif = 1.0 — 2.0. As the correlation U increases, at around 
U ~ 15. OA, the system attains the Kondo regime and the rif assume values close to the 1.0 
as is expected to the Kondo limit of the model. In the inset we represent the limit where the 
/ occupation number becomes exactly rif = 1.0. This point corresponds to the symmetrical 
limit of the Anderson Hamiltonian. 
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FIG. 11: Density of states p f (u) for U = 25. OA, E f = -10. OA and T = 0.001A. 
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FIG. 12: Density of states pf(fJ-) at the chemical potential p as function of the correlation U 
for Ef = —10. OA and T = 0.001 A, for the atomic method and for the Friedel sum rule (FSR). 
The FSR is represented by points over the curve. In the inset we represent the density p c (p) of 
conduction electron states at the chemical potential p, as a function of the correlation U. 

VII. THE DIFFERENT REGIMES OF THE MODEL 

In the set of Figs, fl 1 4ti 1 9 j) we fix the correlation energy value in U = 20. OA and we vary 
the localized level Ef in order to describe the different regimes of the model: Ef = 5. OA 
(empty-dot regime); Ef = 0.0 A (intermediate valence (IV) regime); Ef = —5. OA and 
Ef = —15. OA (Kondo regime), Ef = —20. OA (crossover from the Kondo to the magnetic 
regime) and Ef = —25. OA (magnetic regime). 
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FIG. 13: The total localized occupation number nj as a function of the correlation U in A units. 
In the inset we represent the limit where the / occupation number becomes exactly nj = 1.0. This 
point corresponds to the symmetrical limit of the Anderson hamiltonian. 
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FIG. 14: Density of states for T = 0.001A , U = 20A and E f = 5.0A. This parameter set 
corresponds to the empty dot regime. 

In Fig. [THwe represent the empty dot regime. We have only a tail of the density of states 
below the chemical potential fj, and the total occupation number is very low and is given by 
n f ~ 0.12. 

In Fig. [TB] we represent a typical intermediate valence situation. In this case the density 
of states already exhibit the three peak structures characteristic of the Kondo regime, but 
the Kondo peak is not yet formed; there is a large structure at the chemical potential //, 
that generates a strong charge fluctuation. The total occupation number is characteristic of 
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FIG. 15: Density of states for T = 0.001A , U = 20A and E f = 0.0A. This parameter set 
corresponds to the intermediate valence regime (IV). 
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FIG. 16: Density of states for T = 0.001A , U = 20A and E f = —5. OA. This parameter set 
corresponds to the beginning of the Kondo regime. 

the intermediate valence regime; nf ~ 0.50. 

In Fig. [TBI we represent the beginning of the Kondo regime for Ef ~ —5. OA and nf ~ 0.96; 
the system remain in this regime until Ef ~ — 0.15A and n/ ~ 1.04 as indicated in Fig. [J71 
In this region, the Kondo peak is well defined and is pinned at the chemical potential fi, and 
in Fig. [22] we present a resume of all regimes. The interesting point that should be stressed 
here, is that the Kondo limit, where the localized occupation number is exactly equal the 
unity (jif = 1.0), is attained in the symmetrical limit of the Anderson model as indicated in 
Fig. UM 
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FIG. 17: Density of states for T = 0.001A , U = 20A and E f = -15.0A. This parameter set 
corresponds to the ending of the Kondo regime. 
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FIG. 18: Density of states for T = 0.001A , U = 20A and Ef = -20A. This region is dominated 
by the double occupation band and the / occupation number is greater than one and corresponds 
to the crossover from the Kondo to the magnetic region. 

For Ef < —15. OA the influence of the double occupation band over the Kondo effect 
increases more and the Kondo peak enlarges its width as indicated in the Fig. [THJ In this 
figure Ef ~ U and rif 1.50 and there is no Kondo effect. We call this region the "magnetic 
region" because as Ef becomes more and more negative there is a competition between the 
Kondo state |11 > and the two magnetic states |14 >,|15 > of the atomic solution of the 
SIAM. We will discuss this point in more details in Figs. [2011221 

Finally in Fig. [19] we represent the limit where (U + Ef) « 0; in this particular case 
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FIG. 19: Density of states for T = 0.001A , U = 20A and E f = -25.0A. This limit is com- 
pletelly dominated by the double occupation band and the / occupation number is close to 2.0 
and corresponds the region where Ef — U < 0. 

the double occupation state is almost completely full and rif ~ 1.88. 




Particle number 



FIG. 20: The sixteen energy levels E n r (in units of A) of the AAS, obtained in Appendix [A] 
(cf. Table III), as a function of the particle number for typical parameters of the Kondo regime: 
U = 20A, E f = -15. OA and T = 0.001A. 

It is now convenient to study in more detail the crossover region from the Kondo to 
the magnetic regime. In Fig. [2U] we represent the sixteen energy levels E n , r = Ej (j = 
1, 16) of the AAS (cf. Table III) as a function of the particle number, for the parameters 
corresponding to the ending of the Kondo regime as indicated in the Fig. [T71 U = 20 A, 
Ef = —15. OA and T = 0.001A. The lines joining them are associated to the corresponding 
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FIG. 21: Intensities of the residues associating to the transition "Kondo- magnetic" as a function 
of Ef (in units of A), corresponding to U = 20A, and T = 0.001A. 

poles Ui (cf. Table IV) of the atomic GF in Eq. (1A2j) that are given in Eqs. (1A5HA13]) 

When Ef < 0.15 the double occupation states become available to the system and there is 
a competition between the singlet that originates the Kondo effect represented by the state 
1 11 > in Fig. [20] and a magnetic state represented in the same figure by the two degenerate 
states 1 14 >,|15 >. 

This point deserves a more detailed study. To do so we consider Fig. [2TJ at Ef — —22. OA 
there is a change of ground state from the two-particle Kondo singlet 1 1 1 > to a three particle 
magnetic doublet (|14 >, |15 >) and we call this point Prm ( "Kondo- magnetic" transition). 
In that figure we present the transitions «4, W12, ^io, an d «i4 connected to the singlet |11 >, 
and their intensities are important in the Kondo region but they decrease strongly when Ef 
becomes closer to the point Prm- For values of Ef < Prm the transitions 144,1X11 and ui$ 
associated to the magnetic state grow up. This region becomes relevant when the RKKY 
interaction is present, like in the problem of several interacting impurities or in the lattice 
case. This interaction is associated to the magnetic transitions in the heavy fermion Kondo 
problem, in which case this transition is generally antiferromagnetic. 

One interesting problem in which the above discussion applies is the double quantum 
dots problem (DQD), which has been recently proposed as a possible realization of the spin 
quantum computer [331]. It constitutes the minimal system for studying a lattice of magnetic 
impurities in a tunable environment. In this system the competition between the Kondo 
effect and the RKKY interaction lead to a second-order quantum phase transition (QPT) 
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but if the electrons can tunnel between impurities this QPT is replaced by a crossover 
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FIG. 23: The total localized occupation number nj as a function of the correlation Et. In the 
inset we represent the limit where the / occupation number becomes exactly nt = 1.0. This point 
corresponds to the symmetrical limit of the Anderson Hamiltonian. 

In Fig. [22]we show the density of states at the chemical potential \x vs. Ef for T = 0.001 A, 
U = 20 A. The two vertical lines separate three main regions: intermediate valence to Kondo 
region and the Kondo region to the magnetic region. From this result, it is clear that in 
the case of finite U the Kondo effect only exists in a limited parameter region of the SIAM. 
Fig. [23] shows that the Kondo behavior manifests itself when the total / occupation number 
rif — Ylio n fi° assume values close to one rif — 1.0; according to the inset of the figure, this 
value of rif corresponds to the symmetrical limit of the SIAM as represented in Fig. [101 
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VIII. CONDUCTANCE OF A SIDE-COUPLED QUANTUM DOT FOR THE FI- 
NITE U CASE 



In this Section we apply the atomic method for finite energy correlation U, to study 
the electronic transport through a quantum wire (QW) with a side-coupled quantum dot 
(QD) 23j. The quantum dot can be occupied from zero to two electrons as a function 
of the chemical potential //. This system has been already studied for U — > oo when the 



double occupation is forbidden 
produces interesting resu 
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31] . The finite U case is much more realistic, and 



ts_35[ for the conductance that should be compared with recent 
37J. 



experimental data 

In Fig. [24] we present a pictorial view of the quantum dot, side-coupled to a ballistic 
channel. 

Quantum dot 




Ouantum wire 

FIG. 24: Pictorial view of a quantum dot side-coupled to a conduction channel. 

Electron transport is coherent at low temperature and bias voltage, and a linear-response 
conductance is given by the Landauer-type formula [31 1 
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2e 2 



dn f 



(53) 



h J \ du 

where rip is the Fermi function and S(u) is the transmission probability of an electron with 
energy fkv. This probability is given by 



S(cj) 



r 2 I (ia 1 2 
1 I ^00 I ) 



(54) 



where T corresponds to the coupling strength of the site to the wire. Gq can be cal- 
culated by the Dyson equation, where V = |0)V(1| + |1)V(0| is the hybridization. The 
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dressed Green's functions at the site can be written in terms of the undressed localized 
Green's functions gw at the QD and of the undressed Green's functions goo of the conduction 
electrons: 

G° 00 = 9oo + 9o VG1 + f&VGSo, (55) 

G a w = 9io + 9ioVG^ + g^VG^. (56) 

Solving this system of equations, and taking into account that the non-diagonal bare con- 
duction Green's functions vanish: gf Q = and g§ x = 0, we can write 

/~ta 000 (K7\ 

°°"(l + ^ 2 <l)' [ } 

where 

fz + D + fi 

0oo = 

and M? l (z) is obtained employing the atomic approach Green's functions given by Eq. ([37 
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FIG. 25: Conductance of the side-coupled quantum dot for U = 20. OA and T = 0.001A. 

In Fig. [25] we calculate the conductance corresponding to the parameters employed in 
Sec. I VI I We show the conductance as a function of Ef for U = 20A and at temperature 
T = 0.001A. We can observe four important regions indicated in the graph. The first 
is the empty-dot region. In this situation, the value of the localized level Ef is positive 
and is located above the chemical potential /i; the localized occupation number goes to 
zero {rif — > 0) and the conductance goes to one (G/G — ► 1). In this situation there is 
practically no destructive interference between the dot and the wire. 
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The next regime is the intermediate valence, which is indicated in the graph by "IV 
region", and the localized level Ef is located around the chemical potential fi = 0. In this 
region the value of the / occupation number varies strongly as a consequence of the charge 
fluctuation, the quantum interference effects start to increase and the Kondo peak appears 
when rif — > 1. 

The next region in the graph shows the effect of the Kondo regime in the side-coupled 
QD system. The occupation number is around one (n/ ~ 1) and the conductance goes to 
zero (G/Gq — > 0). The Kondo peak is formed, and there is a perfect destructive quantum 
interference of tt/2 between the electrons that go around the wire and those that "visit" the 
QD and return to the wire. 

In the last region the / occupation is greater than one, and the double occupation dom- 
inates this regime as indicated in Figs. [T8lfT9l This is the "magnetic" region and we can 
observe a competition between the Kondo effect and the magnetic state, as discussed in the 
text of Figs. I20H221 The occupation number goes to two (jif — > 2) as the Ef becomes more 
negative and the conductance goes again to one (G/Gq — > 1). 

IX. CONCLUSIONS 

Employing the atomic approach, we characterized well the formation of the Kondo peak 
for all the parameters of the SIAM by calculating several curves of density of states, varying 
the value of the coulomb repulsion U, as well as the position of the localized impurity level 
E f . 

As a general conclusion we can say that we have developed a general and simple method 
to calculate the low temperature properties of the Anderson impurity model, and we call 
this method the atomic approach, because the starting point of the method is the zero 
conduction bandwidth limit of the Anderson impurity model. The approach is extremely 
simple, satisfies the Friedel sum rule by construction and is valid, at low temperatures, for 
all the relevant range of parameters and for all the coupling regimes of the SIAM, namely 
for the weak, intermediate and strong correlation regimes of the model. 

The atomic approach is a good candidate to describe nanoscopic systems with correlated 
electrons that present Kondo effect, and it produces excellent results for the occupation 
numbers and for the dynamical properties. Due to its simplicity, low computational cost 
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(the computational time consuming to obtain a density of states curve takes few seconds 
or less) the atomic approach could be applied to study nanoscopic correlated systems, like 
quantum dots. 
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APPENDIX A: THE ATOMIC GREEN FUNCTIONS 

In this Appendix we present the details of the calculation of the exact Green's function of 
the Anderson impurity model in the limit of zero conduction bandwidth. In this limit all the 
hoping contributions are eliminated from the Hamiltonian, because we also take Vf^ a = V 
in Eqs. ( TrTHT~2l) . i.e.: a local hybridization. Transforming the conduction electrons to the 
Wannier representation we then have an independent system at each site of the crystal. In 
this limit there is an isolated metal atom at each site, one of them being the Anderson 
impurity, and the system Hamiltonian can be diagonalized exactly. In the Anderson Hamil- 
tonian, there are four possible occupations (0, |, |, Tl) °f the conduction electrons at any 
site, and at the impurity there are also four possible occupations of the local / electrons: 
(0, +|, — ~, d). At the impurity site we then have a Fock space with sixteen states char- 
acterized by \m, a > as shown in Table HIl and the hybridization mixes states with equal 
particle number and with the same spin component z. We should then diagonalize a 16x16 
matrix which presents a block structure that simplifies the calculation, because the greater 
block is 3x3 and we employ Cardano's formula to solve the associated third degree algebraic 
equation. The results of this calculation are presented in Table IIHI 

To obtain the localized atomic Green's functions of the impurity in the zero width limit, 
we use Zubarev's [38[ equation 

• \ = pn V- exp(-/fe w -i, r ) + exp(-/fe ra , r /) 

n,r,r' ' 

x(n-l,r\ X j>a \n,r') (n,r'\ X] a , \n-l,r), (Al) 
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TABLE II: States of the Anderson impurity in the limit 2D = V = 0. The columns indicates the 
states \n, a >, the energies E Q , the number of electrons n and the spin component S z . 



where Q is the thermo dynamical potential and the eigenvalues E n j and eigenvectors \nj > 
correspond to the complete solution of the Hamiltonian. The final result is the following 

16 

G//.«*( w ) = e ^y_^!_ ) (A2) 
where the Ui poles of the Green's functions are given in Table [IV] 

The residues for the localized electrons are given by 

mi = C o S (0) 2 [e-^ + e"^ 2 + f e"^ 4 + \e~^ E % 
m 2 = sin(<J)) 2 [e-^ + e"^ 4 + fe"^ 2 + fe"^ 6 ], 
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Eigenstates \j >= \n,r > 
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TABLE III: Exact solution (A AS) of the Anderson model in the limit of local hybridization and 
conduction band with zero width. The sixteen eigenstates | n,r) have energies E n>r , where n is the 
number of electrons and S z is the spin component, and we abbreviate \ j >= |n, r > and Ej = E n ^ r 
(j = 1, . . . , 16). We use e„ ir = E n>r — nu, and the other parameters are given by 
A = \{e f - e q ) 2 + W 2 } 1 / 2 - A' = [(e f + U - e q ) 2 + AV 2 } 1 / 2 ; 
tg<j> = 2V/{e q -s f + A); tgO = 2V/(e f + U-e q - A'); 

ai = ^ 2+ 4V^[(E l -2e f -U)^+(E,-2e q )-^ ^ = E i -£ t -U a * Ci = E^2T q a ^ 1 = 9 > 10 ' U ! 

6 1 = arccos y=f ! Q = -|[12F 2 + (e, + £/) 2 + (2e f + U) 2 + {2e q ) 2 - (e q + e f )(2e f + U)- (e q + 
e f )(2e q )-(2e f + U)(2e q )]; 

R=l A {-Z[(e q + e f ) 2 {2e f + U) + (E q + e f ) 2 (2e q ) + {2E f + U) 2 (e q + e f ) + (2e f ^ 

e f ) + (2e q ) 2 (2e f + U)} + I2{e q + e / )(2e / + U){2e q ) + 18V 2 [2( £(? + e/ ) - (2e f + U)- (2e q )} + 2[(e q + 

e / ) 3 + (2e / + ^) 3 + (2e g ) 3 ]}. 
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Ui = 


E2 - 


-Ei = 


Eq — 


Ei 


u 2 = 


E 5 - 


-Ei = 


Eq — 


E2 


U3 = 


E\o 


— E2 






Ui = 


En 


— E2 






u 5 = 


Eg - 


- E2 






Uq = 


E10 


-Ei 






u 7 = 


E12 


— Eq = 


- Eiq 


— En 


u s = 


E12 


— Eg 






Ug = 


E12 


— ElQ 






UlO = 


E12 


— En 






Ull = 


En 


— ElQ 






U12 = 


E14 


— En 






U13 = 


Eg - 


-E4 






u u = 


En 


-E 4 






U15 = 


En 


— Eq = 


- Eiq 


— E\2 


uw = 


Eu 


— Eg 







TABLE IV: Poles of the Green's functions, corresponding to all the possible transitions in the 
atomic solution. 



m 3 = ( e -^3 + e~^«')[(a loS m(0)) 2 + (b w cos^)) 2 ], 
m 4 = ( e -^3 + e -^)[(aiism(0)) 2 + (6 llC os(0)) 2 ], 
m 5 = ( e -^3 + e-P E *)[{agsin{<t))f + (6 9 cos(0)) 2 ], 
mQ = (e-P E4 + e-^ El0 )[(ai cos((f))) 2 + (b 10 sin(<f>)) 2 ], 
m 7 = S in{6f[l(e-P E * + e"^ 12 ) + (e"^ 15 + e"^ 16 )], 
m 8 = (e-PE* + e-^ 2 )[(cgCos(6)) 2 + (a 9 sm(0)) 2 ], 
mg = (e-^o + e ~^)[(c w cos(e)) 2 + (a 10 sin(9)) 2 ], 
m 10 = ( e -^ii + e~^)[(ciiCos(6)) 2 + (a u sm(0)) 2 ], 
m u = ( e -^ 10 + e-P El4 )[{ci sin{6)) 2 + (a w cos(6))\ 
m 12 = ( e -^u + e -^ 15 )[( Cll sm(^)) 2 + (ancos(0)) 2 ], 
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m 13 = (e"^ 5 + e- pE9 )[{a Q cos{<j))) 2 + (6 9 sm(0)) 2 ], 
m 14 = (e"^ 5 + e _/3£;il )[(aiicos(0)) 2 + (&nsm(0)) 2 ], 
m 15 = cos(6)) 2 [f (e"^ 8 + e"^ 14 ) + (e"^ 13 + e"^ 16 )] . 
m 16 = (e"^ 9 + e-^)[{a 9 cos{9)) 2 + (c 9 sm(0)) 2 ], 
and for the electrons c we have 



r cc,at r \ _ _ pQ e X P(-^n-l,r) + exp(-/fe ra , r / 



■X 

iuj s + e n _i r — e r 

n,r,r' 

x (n — l,r| |n, r') (n, r'| ct |n — l,r), (A3) 



16 



G oc«t/ w ) = ^ny-_^_ (A4) 

^ - Mj 

i=l 

and the residues for the conduction electrons are given by 



m = sm(0) 2 [e- /3Bl + e"^ 2 + fe"^ 4 + fe"^ 8 ], 
n 2 = cos((f)) 2 [e-^ + + fe"^ 3 + fe"^ 8 ], 

n 3 = (e~^ + e-^)[(a l0 cos(<P)) 2 + (c lo sm(0)) 2 ], 
n 4 = ( e - pE2 + e" /3£ll )[(ancos(0)) 2 + (c u sm(</>)) 2 ], 
n 5 = (e"^ 2 + e- /3£9 )[(a 9 cos((/))) 2 + (c 9 sm(0)) 2 ], 
n 6 = (e - ^ 4 + e-^ lo )[(a lo sm(0)) 2 + (c lo cos(0)) 2 ], 
n 7 = cos(0) 2 [e-^« + e"^ 14 + fe"^ 8 + fe"^ 13 ], 
n s = ( e -PE<> + e -^ 2 )[(a 9 cos(^)) 2 + (6 9 sm(0)) 2 ], 
n 9 = ( e -^io + e -^)[(a w cos(9)) 2 + (& lo sin(0)) 2 ], 
n 10 = (e"^ 11 + e- /3Bl2 )[(a 11 cos(^)) 2 + (b u sin(9)) 2 ], 
n u = (e-^ 10 + e- /3£l4 )[(ai sm(^)) 2 + {b l0 cos{6)) 2 }, 
n 12 = (e"^ 11 + e-P E ^)[{b u cos{6)) 2 + {a u sin{6)) 2 ]. 
n 13 = (e^ 4 + e _/3&J )[(c 9 cos(0)) 2 + (a 9 sm(0)) 2 ], 
n 14 = (e - ^ 4 + e _/3£;il )[(ciicos(0)) 2 + (a n sm(0)) 2 ], 
n 1B = sin(e) 2 [e-^ + e"^ + fe"^ 14 + fe"^ 6 ], 
n 16 = (e-^ 9 + e- /3El4 )[(6 9 cos(^)) 2 + (a 9 sm(0)) 2 ], 



Finally following the same definition for the cumulants in Eq. (123|) we can write all the 
atomic Green's functions employed in the calculation of the atomic approach. 
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GffM _ 




Qff,at 



922 924 
942 944, 



911 = 6 { [ Sm U + + 2 l u + E 2 -E & ) 

( e -t>Ei + e ~PE 2 3 e -/3£ 4 + e -/3£ 6 X "I 

+ cos 6 1 

Y V + E x - E 2 2itu + E i - E 6 J _ 

4^ iu; + 



i=9 



i _|_ g— /3Ei4 

icj + Ei 



E- (a - Sin0)2+ zc + E 5 -^ (a - 



'14 



(c* sinfl) 2 + 



— (Oj cos0) 2 

e ~PEi _|_ e -/3^i2 



- - - v qcos6>) 2 



/ e -/3£i5 _)_ e -/3Bi6 3 g-/3S 8 _|_ e -/3Ei2 

sm 2 6» 1 — 

\ iu + Eva - E 16 2 iu + E 8 - E12 

/ e -f3E 13 I e -/3E 16 3 g-/3£ 8 _i_ e -/3Si4 \ " 

+ C ° S2 6 ET- + o — — ~ 

11 - 



p -/3£ 4 i p-PE, -flE z , 

~ / 7 • / \ 2 ~ / 7 / \ 2 

e -i3E 12 „-0Ei , „-/3Bi4 1 ^ 



2iu + E s - E lA 



+ 



e -f3Ei _|_ e -/3Ei2 g-/3£i _|_ e -/3Ei4 

— (di sin#) 2 H — (cii cos9) 2 

iu + Ei- E 12 K ' iu + Ei - E 14 K ' _ 



+ 



n 

E 

i=9 
-0E, 



+ e 



-(SEi 



-/3E 3 



+ e 



iu + E 5 - Ei iu + E 3 - Ei 



(dfii sin <fi cos 0) 



+ e 



-/3Bl4 



iu + Ei- E 



14 



12 



fa,c,- sin 6 cos i 



531 — 9l3, 



sin 2 (j) 



-PEi 



-PE 5 



+ cos 4> 



iu + Ei — E 5 

PEl _|_ g-/3£ 3 



+ 



3e 



2 icj + E-x - E 7 



+ 



3 e ~PE 5 + g-/3S 7 



iu + E ± - E 3 2iu + E 5 - E 7 



li 



+ e" 



+ 



i=9 

e -/3£i + g 



■(oj sin 0) 2 + 



+ e- 



(a, COS0)^ 



+ - £"15 



-(Qsin^)^ + - 



iu + E 4 — E, t 

(ci cos 9) 2 



e —PEi _|_ g— (3Ei3 



iu + Ei - E 
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944 



+ cos 2 6 
11 



sin 9 



e -f3E 14 _|_ e -0E 16 



iu + E\a — E 



+ 



3 e -PE a + e -(3E l2 



+ e 



ii 



n6 



2iu + E 6 - E- 



12 



+ 



E 



i=9 

4 



icu + Exi - E m 
{hi sin i 



3 e -/3E a _|_ g -/3£;i4 



2 icj + E G - E 



11 



zu; + E4 — E, 



+ 



e -/3E 2 _|_ e -(3Ei 

iuj + E 2 - Ei 



(bi cos 4>y 





zo; + - £13 



a, sm 1 



+ 



+ e 



iuj + Ei — E15 



di cos ( 



(All) 



924 



+ 



11 

E 

i=9 
-(3E 



-/9J3 4 



iu + E 2 — Ei 



+ e 



g— /3-Ei _|_ g— /3_Ei5 



IU + Ei~ ^13 



icj + £i - E 15 



ia; + £?4 — E{ 

aid sin # cos f 



(a,;&,: sin cos < 



(A12) 
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